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The scaling of the tails of the probability of a system to percolate only in the 
horizontal direction tt^s was investigated numerically for correlated site-bond perco- 
lation model for q = 1,2,3,4. We have to demonstrate that the tails of the crossing 
probability far from the critical point p c have shape 7r^ s (p) — D exp(cL[p — p c ] u ) 
where v is the correlation length index, p = 1 — exp(— /?) is the probability 
of a bond to be closed. At criticality we observe crossover to another scaling 
TThs(p) — ^4exp (— b{L\p — p^Y). Here z is a scaling index describing the central 
part of the crossing probability. 



INTRODUCTION 



The g-state Potts model can be represented as the correlated site-bond percolation in 
terms of Fortuin-Kasteleyn clusters lj. At the critical point of the second order phase 
transition, the infinite cluster is formed. This cluster crosses the system connecting the 
opposite sides of the square lattice. In the last decade the study of the shape of the crossing 
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as well as numerically 



the distribution function of the percolation thresholds is 



Gaussian function. Following the number of works 
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the tails of the 



distribution function are not Gaussian ones. The authors of the recent work Ref. 
still uncertain to distinguish a stretched exponential behavior from a Gaussian. 
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The aim of this paper is to investigate the shape of the probability of a system to percolate 
only in horizontal direction 717^. We perform numerical simulation of correlated site-bond 
percolation model for q = 1,2,3,4 (the percolation model q — 1, the Ising model q = 2 
and the Potts model q = 3,4) for lattice sizes L = 32,48,64,80, 128. The scaling formulas 
for a body of the crossing probability at criticality and for tails of the crossing probability 
were obtained. The final result for the representative case q = 2, L = 128 is immediately 
presented in Fig. [TJi). Details of fitting procedure are described in Section El In this figure 
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FIG. 1: a) The absolute value of the logarithm of the crossing probability — ln(7r^ s ) for the Ising 
model as a function of the absolute value of the scaling variable r = L(p — p c ) u on the lattice 
L = 128; b) The magnetic susceptibility x(P) an d the crossing probability vr/j(/3) as a functions of 
the inverse temperature (5 for the Ising model, L = 128. 
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we plot (by crosses) the numerical data for the absolute value of the logarithm of crossing 
probability 717^ for the Ising model (q = 2) on the lattice L = 128 as a function of the 
absolute value of the scaling variable r = L(p — p c ) u - Here p = 1 — exp(— 0) is a probability 
of a bond to be closed, (3 is the inverse temperature, v is the correlation length scaling index. 
The critical point in the p scale for the g-state Potts model is p c = ^/q/(l + ■s/q) see Ref. 2^| 
and we get p c (q = 2) ~ 0.58578 

We can see from Fig.[IJi) that the function 7r^ s (r) consists of two parts: the body |t| < r 
and the tails |r| > r . The negative logarithm of the body of the crossing probability as a 
function of r is well described by function — \n(A) + b\r\ z (solid line on the Fig. UK)). Here z 
is some scaling index. The value of the crossing probability at the critical point TThsiPc) — A 
may be computed (at least for percolation) by conformal field methods 4]. The negative 
logarithm of the tails of the crossing probability have shape — ln(D) + cr (dashed line in 
the Fig. [T^)). This line is tangent to the body at the point r . This point is marked by the 
horizontal line. Let us note that in Fig. QJi) we plot two branches of the crossing probability 
(for r > and r < 0). The coincidence of this two branches indicate the remarkable 
symmetry of the crossing probability with respect to the variable r. 

In Fig. QJd) we plot the crossing probability by crosses (bottom) and the magnetic sus- 
ceptibility by triangles (top) as a functions of the inverse temperature (3 with logarithm 
scale for the ordinate axis. In Fig. [TJa) we indicate the position of crossover region of tt^ 
by horizontal solid line on a level D exp(— ctq). For the magnetic susceptibility we mark the 
region with critical behavior 

x(0)~((0-&)/&r r (!) 

by horizontal dashed lines. We see from Fig. QJd) that tails of the crossing probability directly 
correspond to the critical region of the magnetic susceptibility. In this critical region the 
correlation length £ is smaller than the sample size £ < L. As the temperature approaches to 
the critical point, the correlation length reaches the sample size. At that point the magnetic 
susceptibility on the finite lattice deviates from the critical behavior Eq. ([!} and becomes 
smooth - see the region over the top dashed horizontal line in the Fig. [1J)). At the same 
point the crossing probability crosses over from tails to body - the region over the solid 
horizontal line in Fig. Qt) (and the region under the horizontal line in Fig. QJi))- At the 
critical point (3 C = — ln(l — p c ) ~ 0.881373... both the magnetic susceptibility and the 
crossing probability reach a maximum. 
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The detailed description of the fitting procedure as well as numerical data for q = 1, 2, 3, 4 
are presented below. The main numerical result of this paper is the proving of the formula 
D exp(b[p—p c ]L 1 ') for the tails of the crossing probability. Therefore, we pay special attention 
to fitting procedures. 

The paper is organized as follows: In the second section, we describe details of the 
numerical simulation. In the third section, the method for determining the pseudocritical 
point p c {L, q) on the finite lattice is described. We use p c (L, q) to perform the approximation 
of the tails. In section IIVI we approximate the double logarithm of the crossing probability 
ln(— ln(7Ths{p', L, q))) tails as a function of the logarithm of deviation from the critical point 
ln(p — p c (L, q)) by the linear function c(L, q) +x(L, q) \ia(p — p c (L, q)). We get x(L, q) ~ u(q) 
for this approximation procedure. In section E] we describe new fitting procedure using the 
scaling variable r = L(p — p c ) u . Results of approximation are discussed in Section IVTl 



II. DETAILS OF NUMERICAL SIMULATION 



We perform the massive Monte-Carlo simulation on the square lattice of size L to obtain 
the high-precision data for 717^. We use the dual lattice shown in Fig. |21 On such lattice 
the critical point of the bond percolation (q = 1) is exactly equal 1/2 and is not dependent 
on the lattice size 2^|. To produce the pseudorandom numbers we use the R9689 random 
number generator with four taps 25]. We close each bond with a probability p and leave 
it open with a probability 1 — p. Then we split the lattice into clusters of connected sites 
by using the Hoshen-Kopelman algorithm [26j • After that we check the percolation through 
this configuration. We average the crossing probability over 10 7 random bond configura- 
tions. To investigate the Potts model q = 2, 3, 4 we use the Wolff j^f cluster algorithm to 
generate a sequence of thermally equilibrated spin configurations. For each particular in- 
verse temperature (3 = ^ we flip 20000 Wolff clusters to equilibrate the system. For the spin 
model on the finite lattice the deviation of the pseudocritical point from the position of the 
critical point on the infinite lattice is smaller for Periodic Boundary Condition (PBC) rather 



than the Open Boundary Condition (OBC) 28|, |29|. For this reason we use the PBC for the 



Wolff algorithm. The Monte-Carlo algorithm generates spin configurations on a torus. For 
a generated spin configuration we create a configuration of bonds. Each bond between sites 
with equal spin variable a is closed with the probability p — 1 — exp(— (3) and is open with 
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FIG. 2: The dual lattice with the horizontally spanning cluster. 



probability l—p = exp(— (3). Bonds between sites with different values of a are always open 
in accordance with Fortuin-Kasteleyn rule [l]. Then we split the particular spin and bond 
configurations into different clusters. Here we use OBC It means that for each generated 
configuration we cut the torus and check the crossing on the square with open boundary 
conditions. We fix the OBC for crossing only in horizontal direction 717^ because it implies 
the vertical crossing is absent and the top and bottom rows must be disjointed. But we take 
into consideration the additional raw and column of bonds, as shown in Fig. El We check 
the percolation through an obtained cluster configuration, generate new spin configuration 
and so on. We average the crossing probability over 10 7 configurations for each value of the 
inverse temperature (3. So the resolution of our computations is about 10~ 7 . By this way 
we perform numerical simulation and get a set of data for 717^ (p; L,q) for the lattice sizes 
L = 32,48,64,80,128 and q = 2,3,4. The formal definition of irhs(j>',L,q) as a sum over 
different cluster configurations is described in j^ . 

For the Potts model we use the dual lattice as we do for the percolation. It means, 
that we take into account additional bonds attached to the bottom row of spins. In the 
same way we take into account additional bonds attached to the right column of the spins. 
On the lattice with PBC these bonds have to connect the right and the left columns. We 
cut the torus (because we use OBC for crossing probability) but we keep these additional 
bonds and take into account the checking of the crossing. In Fig. 121 contact points are 
shown by arrows. The left contact points are attached to the left column of sites. The right 
contact points are attached to additional bonds. In Fig. |21 the bond configuration with the 
horizontally spanning cluster is shown. Then we check the percolation through the obtained 
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cluster configuration. After that we flip three Wolff clusters, check spanning for a new spin 
configuration and so on. 

III. DETERMINATION OF THE PSEUDOCRITICAL POINT ON THE FINITE 

LATTICE 

We investigate the crossing probability as a function of deviation from the critical point. 
Therefore, we perform the preliminary approximations to obtain the critical points for the 
finite samples. Namely, we: 

• obtain the pseudocritical point and the shape of the central part of the crossing prob- 
ability, 

• determine the shape of the tails of the crossing probability 

• combine together the information for body and tails, and reconstruct the total shape 
of the crossing probability. 

We need to recall that we consider the crossing probability as a function of the variable 
p — 1 — exp(— (3) (probability of a bond to be closed). It is easy to understand that we must 
take the pseudocritical point on the finite lattice p c (L, q) as the reference point. The crossing 
probability is a symmetric function of the variable Ap = (p— p c (L, q)). This fact implies that 
the high temperature tail p c — Ap (Ap > 0) and the low temperature tail p c + Ap coincide 
KhsiPc ~ Ap) = TThs(p c + Ap). For the bond percolation on the dual lattice the position of 
the percolation point does not depend on the lattice size p c (L, q = 1) = 0.5 j^J. 

For the determination of vJL, q) for the Ising and the Potts model we use the following 
procedure. We can assume [2QJ that in the region —6 < ln(7T/j S (p; L, q)) the fitting formula 
is true 

F(p; L, q) = A(q, L) exp(- [B(L, q)(p - p c (L, q))f L ' g) ). (2) 

Therefore we fit the logarithm of the crossing probability ln(7r^ s ) by the function f\(p; L, q) = 
ln(F (p; L , q)) namely 

A(p; L, q) = a(L, q)) - [B(L, q)(p - p c (L, q))f L ^ , (3) 

where a(L, q) = ln(v4(L, q)). We plot the data for ln(7r^ s (p; L, q)) as a function of p in Fig.^K), 
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FIG. 3: Approximation of the logarithm of the crossing probabilities In (717^) as a function of p by 
the function f\{p;L,q) = a(L,q) — B(L,q)(p — p c {L,q))~^ L ' q ^> for: a) percolation q = 1, b) Ising 
model q = 2, c) Potts model q = 3, d) Potts model q = 4. 

Fig. Eb), Fig. Efc), Fig. Eli) for g = 1, 2, 3, 4 respectively. The errorbars in these figures are 
about the symbol size. It seems, that behavior of \n(irh s {p', L, q)) near p c is parabolic. Results 
of the approximation are plotted in the same figures by lines. We shall see that there is 
a good agreement between the numerical data and the results of the approximation. But 
we see deviation at the point p c especially for q = 3,4. In the vicinity of p c real graphs 
are more smooth than fitting functions. Finally for each pair of numbers (L, q) we obtain 
four fitting parameters A(L,q), B(L,q), p c (L,q), ((L,q). Here A(L, q) defines the crossing 
probability in the critical point, B(L, q) is a scaling variable, p c {L, q) is the position of the 
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pseudocritical point on the lattice L, and ((L, q) is a scaling index. 

In Table H] we collect data for the logarithm of critical amplitude a(L,q) = \n(A(L, q)). 
The fitting parameter a(L, q) defines a vertical shift of graphs in Fig. EH)-IHtO from the zero 
level. In Table HT1 the data for the scaling variable B(L,q) are shown. We approximate the 



TABLE I: Results of the approximation for the fitting parameter a(L,q) = ln(A(L,q)). 



L 


q = l 


q = 2 


q = 3 


q = 4 


32 


-1.708(4) 


-1.773(14) 


-1.881(23) 


-2.022(27) 


48 


-1.710(4) 


-1.773(16) 


-1.886(27) 


-2.043(30) 


64 


-1.711(4) 


-1.776(16) 


-1.894(27) 


-2.055(32) 


80 


-1.711(5) 


-1.777(17) 


-1.897(27) 


-2.078(29) 


128 


-1.705(8) 


-1.777(18) 


-1.898(29) 


-2.089(31) 



data for B(L,q) by the function b*(L,q)L u ^ L ' g \ In Fig. Hk) the results of approximation 
are plotted by lines. Values of the fitting parameters b*(L, q), u(L, b) as well as the inverse 
correlation length index \ are placed in the three last rows of Table HI1 We can see for 
q = 1, 2, 3 that u(L, q) 



K<z) 



. It can be assumed that 



B(L,q) = b*(L,q)L^n. 



(4) 



In the case q = 4 the scaling index u(q) is not equal 



31 



, 1 . Many critical quantities in the 
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These logarithmic correc- 
1.5 and numerical approxima- 



Potts model q = 4 exhibit logarithmic corrections 
tions explain the difference between analytical value ^t~^ 
tion for the scaling index u = 1.372(8). The position of the critical point as a function of 
the lattice size L is shown in Fig. 0Jd). On the dual lattice the position of the critical point 
for percolation is equal | and does not depend on the lattice size. For our computation for 
all points p c (L,q) the deviations from 0.5 are less than 0.0001. This deviation corresponds 
to numerical inaccuracy of our Monte-Carlo simulation. The data p c (L,q) for q = 2,3,4 
were approximated by the power function of the lattice size p c (L,q) ~ p c (<?) + dp(q)L v( - q K 
Results of approximation are placed in Table ILTTl and are also plotted by lines in Fig. ^p). 
We shall see that our fitting procedure determines the critical point with accuracy up to 
four digits after the decimal point. Thus, we shall use obtained values p c {L, q) for follow- 
ing approximation of tails of the crossing probability. In Table IIVI we place the results of 
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FIG. 4: a) The fitting parameter B(L,q) as a function of L. Results of the approximation by 
functions b*(L,q)L u ( L ' q ^ (see Table ITT)l are plotted by lines; b) Position of the pseudocritical point 
on the finite lattice p c (L, q) as a function of the lattice size L for q = 1, 2, 3, 4 and approximation 
by the pow law p c (L, q) = p c (oo, q) + dp(q)L v ( q > for q = 2, 3, 4. Line 0.5 is also added. 



the approximation for the index ((L,q) describing the curvature of the central part of the 
crossing probability. 
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TABLE II: Data for B{L,q) and fitting parameters b*(L,q), u(L,q). 



L 


q=l 


q = 2 


q = 3 


q = 4 


32 


24.8(1) 


60(1) 


116(3) 


195(6) 


32 


24.4(1) 


57(1) 


105(2) 


173(4) 


48 


33.0(1) 


86(1) 


171(4) 


301(8) 


64 


40.8(1) 


114(1) 


240(5) 


449(12) 


80 


48.2(1) 


143(2) 


313(7) 


599(15) 


128 


68.9(4) 


229(3) 


553(13) 


1163(31) 


b*(L,q) 


1.84(2) 


1.767(9) 


1.63(2) 


1.49(5) 


u(L,q) 


0.746(3) 


1.002(1) 


1.198(4) 


1.372(8) 




0.75 


1 


1.2 


1.5 



TABLE III: Results of the approximation of p c (L,q) by the power low p c {q) + dp(q)L<l\ The 
analytical values pP recise (q) = ^/q/(^/q+ 1) are added for comparison. 



q 


dp(q) 


v(q) 


Pc(q) 


pP recise (q) 


2 


-0.160(5) 


-1.25(1) 


0.58573(2) 


0.585786... 


3 


-0.22(2) 


-1.33(2) 


0.63394(2) 


0.633975.. 


4 


-0.26(2) 


-1.42(4) 


0.66666(1) 


0.6666666.. 



IV. THE SHAPE OF THE THE CROSSING PROBABILITY n HS TAILS: 

DIRECT APPROACH. 

Let us check the shape of the crossing probability tails. The double logarithm of the 
crossing probability ln(— ln(7r^ s (L, q))) are plotted as functions of the variable t = ln([p — 
p c (L,q)\) in Fig. E^), Fig. Eb), Fig. Et), Fig. Eli) for q = 1,2,3,4 respectively. Points for 
high-temperature and low-temperature tails lay on the same lines. This yields that we shall 
prove one more time our fitting procedure. 

We expect that the crossing probability tails are described by the formula 

n hs (p; L, q) = D(L, q) exp (-C(L, q)(p - p c (L, q))^) . (5) 
In the interval ln(6) < ln(— In (717^ (L, q))) < ln(12) the tails of the crossing probability look 
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TABLE IV: Results of the approximation for 



L 


q = l 


q = 2 


q = 3 


q = 4 


32 


1.887(6) 


1.52(2) 


1.38(3) 


1.28(3) 


48 


1.883(7) 


1.52(2) 


1.37(3) 


1.27(3) 


64 


1.882(6) 


1.52(2) 


1.38(3) 


1.28(3) 


80 


1.885(7) 


1.52(2) 


1.38(3) 


1.30(3) 


128 


1.86(1) 


1.52(2) 


1.37(3) 


1.28(3) 



like straight lines. 

The absence of prefactor before the exponent in Eq. (0) was argued in Ref . |l2j for the case 
of wrapping in the horizontal direction in terms of transfer-matrix. We can not directly apply 
the same arguments in our case, because we consider another function - the probability of 
crossing only in the horizontal direction. The presence of prefactor must cause the deviation 
from the linear dependence in Fig. El in accordance with Eq. (jUJ) 

In [- \n(ir hs )} = In [—d(L, q) + C(L, q) exp(tx(L, q))} ~ 

(oj 

^HC(L,q))+x(L,q)t- c{ J^ iLq)) . 
Here d(L,q) = ln(D(L,q)) and t = ln(\p — p c (L,q)\) is a logarithm of deviation from the 
critical point. But deviation due to prefactor D(L,q) ^ 1 exponentially decreases as t 
growth so it is possible to avoid the deviation from the linear dependence by appropriate 
choice of an interval of approximation. In the region of approximation points in Fig. [S]lie on 
lines with good accuracy. We obtain the set of p c {L, q) for L = 32, . . . , 128 and q = 2, 3, 4 in 
Section ITTT1 Using this data let us approximate the double logarithm of tails ln(— ln(7r/j S )) 
as a function of t by formula (J7|) 

f 2 (t;L,q) = c(L,q)+x(L,q)t. (7) 

Combining (J3J) and (J7J), we obtain c(L, q) = ln(C(L, q)). The resolution of our computations 
10~ 7 is about 16 units in (—In) scale and is about 2.7 units in (ln(— In)) scale. We use the 
points in the interval ln(6) < ln(- \n(n hs )) < ln(12), ln(6) ~ 1.79, ln(12) ~ 2.48 for this 
approximation. This interval is indicated in Fig. Et)-Fig. EJi) by horizontal lines. In this 
figures the slope lines represent results of approximation for x(L, q). 

In the region ln(— ln(7T^ s )) < ln(6) the crossing probability obeys another scaling formula. 
In the region ln(— ln(7r^ s )) > ln(12) the numerical inaccuracy becomes large. There are 
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FIG. 5: Approximation of the double logarithm of the crossing probabilities ln(— ln(7T/j S (p; L, q))) 
by the function /2(i; A <?) = c(L, q) + x(L, q)t of the logarithm of the distance to the critical point 
t = ln(|p — p c (L, q)\) for: a) percolation q = 1, b) Ising model g = 2, c) Potts model g = 3, d) Potts 
model q = 4. The range of the approximation is shown by horizontal lines. 

only a few numbers of "hits" in horizontally spanning clusters for ln(— ln(7T/j S )) > ln(12). 
Results of the numerical approximation for x(L, q) and c(L, q) are presented in Table IVI 
and Table IVT1 respectively. As we can see x(L,q) ~ v(q). Exception is the case q — 4. 
For q = 4 we obtain x(L,q = 4) ~ 3/4 instead of u(q = 4) = 2/3. This deviations can 
be explained by the logarithmic corrections. Some deviations exceeding the approximation 
errors can be explained by the choice of the approximation region. As it can be seen from 
Fig. EK)-Fig. decreasing of the bottom approximation bound ln(6) causes decreasing of 
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TABLE V: The slope of lines: results of the approximation for x(L,q), the analytical values v(q) 
are presented for comparison. 



L 


q=l 


q = 2 


q = 3 


q = A 


V 


4/3 ~ 1.3333 


1 


5/6 ~ 0.83333 


2/3 ~ 0.6666 


32 


1.366(11) 


1.026(23) 


0.87(2) 


0.761(21) 


48 


1.375(10) 


1.030(12) 


0.868(12) 


0.772(14) 


64 


1.376(12) 


1.035(8) 


0.877(8) 


0.767(10) 


80 


1.32(3) 


1.032(10) 


0.872(13) 


0.752(14) 


128 


1.31(4) 


1.036(5) 


0.874(6) 


0.759(4) 



the slop of the approximation line. From all said above we can conclude that the tails of 
the crossing probability behaves like exp(— (p — p c ) u ). We fit the parameter c(L,q) by the 

TABLE VI: Shift of the lines: results of the approximation ln(C(L, q)) = c(L, q) = d[ + c\ ln(L). 



L 


q=l 


q = 2 


q = 3 


q = A 


32 


5.09(2) 


4.96(6) 


4.940(8) 


4.88(8) 


48 


5.52(2) 


5.39(4) 


5.35(5) 


5.35(6) 


64 


5.82(3) 


5.70(3) 


5.68(3) 


5.63(4) 


80 


5.9(8) 


5.92(4) 


5.90(5) 


5.79(7) 


128 


6.33(12) 


6.42(2) 


6.39(3) 


6.32(2) 



expression c\ + c 2 ln(L) and show results in Table IVTTl In Fig. El the fitting parameter c(L, q) 
is plotted as a function of L for q = 1,2,3,4. Results of approximation c\ + Cgh^L) for 
q = 1,2,3,4 are added as well. We may assume that the parameter c* 2 for q — 1,2,3,4 is 
equal 1 within accuracy of the approximation. In accordance with scaling theory for the 
system of size L the deviation from the critical point is described by expression L(p — p c ) v . 
Therefore we may assume that 

n hs (p; L, q) ~ D(L, q) exp (-c(g)L [p - p c {L, q)f q) ) . (8) 

In Fig. QJi) we show the crossover from the parabolic like functional dependence in the 
vicinity of the critical point to the tails. Namely, we plot the data — ln(7T/j S ) for q = 2, 
L = 128 as a function of the variable r = L(p — p c (L,q)) u . We shall call the function 
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3.4 3.6 3.8 



4 4.2 4.4 4.6 4.8 
ln(L) 



FIG. 6: The fitting parameter c{L,q) = \n(C(L,q)) (the vertical shift of the approximation lines 
for tails) as a function of m(L) for q = 1,2,3,4. Results of approximation c\{q) + c| ((z) ln(L) are 
shown by lines. 



TABLE VII: Results of the approximation of c(L, q) by the function c{{q) + c^iq) ln(L). 



q 


1 


2 


3 


4 


c*i 


1.7(2) 


1.32(12) 


1.30(6) 


1.38(1) 


4 


0.98(6) 


1.05(2) 


1.05(2) 


1.02(2) 



under consideration in the vicinity of the critical point a "body" of the crossing probability. 
Procedure of approximation for this figure is described below. 



V. NEW FITTING PROCEDURE 



From the above said we can make the following conclusions. There are two scaling regions: 
the first one is in the vicinity of the critical point, and the second is the tails. Analyzing the 
Eqs. (J2J), (Jl|) and (jEJ) we can conclude that the distance from the critical point (p — p c (L, q)) 
and the lattice size L occur in a formula only as combination L/£ ~ L(p — p c (L, q)) u ^ in 
accordance with scaling theory. Let us introduce a scaling variable 

r = L(p-p c (L,q)y^. (9) 

If we plot the negative logarithm of the crossing probability — ln^s) as a function of 
the scaling variable r then we expect the power dependence in the vicinity of zero and 
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linear dependence for tails, as we can see in Fig. ^i). We can use the scaling variable 
t = L(p — p c (L, q)) v to fit the crossing probability taking into account the finite size 
scaling. Let us describe a new fitting procedure. 

We may assume that we obtain the position of the critical point on the finite lattice 
p c {L, q) as a result of the previous fit. Thus, we can use only three free fitting parameters 
and fix the value p c (L,q). Substituting Eq. (jHJ) for r in Eq. (jSJ) and Eq. (JSJ), we obtain the 
fitting formulas for the body and the tails of the crossing probability 

- lnK,) ~ / 3 (r; L, q) = -a{L, q) + b(L, q)r< L ^ , (10) 

- ln(7r fts ) ~ / 4 (t; L, q) = -d(L, q) + c(L, q)r. (11) 

Here r = L(p - p c (L, q))^ q) . If our assumptions are correct then fitting parameters a, d, c, z 
do not depend on the lattice size L. 

The numerical value of the variable p c {L, q) is used for computation of r in accordance 
with Eq. @. Comparing new fitting procedures Eqs. (|l()jl. (fTTjl and previous fitting proce- 
dures Eqs. (J3J), (0) we can obtain relations between new fitting parameters b(L,q), z(L,q), 
c(L,q) and parameters b*(L,q), ((L,q), C(L,q) 

b(L,q) = b*(L,q)^\ (12) 

((L,q) = z(L,q)u(L,q), (13) 

C(L, q) = c(L, q)L, c(L, q) = ln(C(L, q)). (14) 

The special case is q = 4. As we can observe from the results of the approximations in 
Table 10] and Table 13 the finite size scaling of the crossing probability for q = 4 does 
not obey the scaling low L(p — p c Y^ q=/ ^ (probably to logarithmic corrections to scaling). 
Therefore, instead of the analytical value u(q = 4) = 2/3 we use the numerical approximation 
x(q = 4) = 3/4 from Table IVl We can see later, that this substitution allows us to obtain 
correct results for scaling. 

We perform the fitting procedure in accordance with the formula proposed above. For 
the fit of the body of the crossing probability Eq. (JTHJ) the scaling region 1.6 < — hi(Tih s ) < 6 
was used. For the fit of the tails of the crossing probability Eq. (jllj) we use the scaling region 
6 < — \n(7i hs ) < 12. We place results of the fitting procedure for q = 1,2,3,4 in Fig. Uh), 
Fig. |7|d), Fig. [Tt), Fig. |7|i) respectively. In these figures the fitting regions are denoted by 
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L=32 
L=48 
L=64 
L=80 

L=128 
f 3 (x;32;l) 
f 3 (X;48;l) 
f 3 (x;64;l) 
f 3 (x;80;l) 
f 3 (x;128;l) 
f 4 (x;32;l) 
f 4 (x;48;l) 
f 4 (x;64;l) 
f 4 (x;80;l) 




-32 

L=48 
L=64 
L=80 

L=128 
f 3 (x;32;3) 
f 3 (x;48;3) 
f 3 (x;64;3) 
f 3 (x;80;3) 
f 3 (x;128;3) 
f 4 (x;32;3) 
f 4 (x;48;3) 
f 4 (x;64;3) 
f 4 (x;80;3) 
f 4 (x;128;3) 




(b) 



L=32 
L=48 
L=64 
L=80 

L=128 
(x;32;2) 
x;48;2) 
,(x;64;2) 
f 3 (x;80;2) 
f 3 (x;128;2) 
f 4 (x;32;2) 
f 4 (x;48;2) 
f 4 (x;64;2) 
f 4 (x;80;2) 
f 4 (x;128;2) 




L=32 
L=48 
L=64 
L=80 

L=128 
f 3 (x;32;4) 
f 3 (x;48;4) 
f 3 (x;64;4) 
f 3 (x;80;4) 
f,(x;128;4) 
f 4 (x;32;4) 
f 4 (x;48;4) 
f 4 (x;64;4) 
f 4 (x;80;4) 
f 4 (x;128;4) 



FIG. 7: The negative logarithm of the crossing probability — ln(7r/, s (p; L, q)) as a function of the 
the scaling variable r = L(p — p c (L,q)) u ^ for: a) percolation q = 1, b) Ising model q = 2, c) 
Potts model q = 3, d) Potts model q = 4. Results of the approximation of the body of the crossing 
probability by the function fs(T;L,q) = —a(L,q) + b(L,q)r z and tails of the crossing probability 
by the function /4(r;L,g) = —d(L,q) + c(L,q)r are added. The ranges of approximation regions 
are shown by horizontal lines. 



horizontal lines. As we expect, if we plot the data as functions of the scaling variable r then 
the finite size dependence eliminates and points for different values of L lie on the same 
curves. Results of the approximation for the fitting parameters b(L,q), z(L,q), c(L,q) and 
d(L,q) are collected in Table IVTTH Table ITX1 Table IX| Table IXI1 respectively. 

The data for b(L, q) in Table IVIIII demonstrates, that this fitting parameter does not 
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depend on the lattice size L. All dependence on L is enclosed in the expression for r. Hence 
we can omit the variable L in the round brackets and consider b(q) only as a function of q. 
This fact proved the choice of the fitting procedure. The data in Table HX] presents some 

TABLE VIII: Results of the approximation for the fitting parameter b(L,q). 



L 


9=1 


9 = 2 


9 = 3 


9 = 4 


32 


3.063(3) 


2.38(2) 


1.91(2) 


1.87(3) 


48 


3.032(3) 


2.40(2) 


1.89(2) 


1.95(2) 


64 


3.018(5) 


2.38(2) 


1.90(2) 


1.97(3) 


80 


3.027(5) 


2.39(2) 


1.91(3) 


2.02(3) 


128 


3.049(6) 


2.38(2) 


1.90(3) 


2.05(3) 



power z(q). This power describes the behavior of the crossing probability in the vicinity of 
the critical point as a function of the scaling variable r and does not depend on the lattice 
size. The data for the fitting parameter c(L, q) are placed in Table [Xj It seems, that this 

TABLE IX: Results of the approximation for the fitting parameter z(L,q). 



L 


9 = 1 


9 = 2 


9 = 3 


9 = 4 


32 


1.432(3) 


1.61(2) 


1.75(3) 


1.82(3) 


48 


1.426(4) 


1.59(2) 


1.77(3) 


1.78(3) 


64 


1.413(5) 


1.59(2) 


1.76(4) 


1.78(4) 


80 


1.432(4) 


1.60(3) 


1.73(4) 


1.73(3) 


128 


1.432(6) 


1.61(3) 


1.74(4) 


1.81(5) 



parameter does not depend on L. The results for d(L, q) is placed in Table IXfl First of all, 
this parameter does not depend on the lattice size within accuracy of the approximation. 
The absolute value of parameter d(q) is relatively small so the prefactor D(q) is about 1. 
For percolation the order of d(q = 1) is less than three value of numerical inaccuracy, thus 
it is possible d(q = 1) = 0. 
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TABLE X: Results of the approximation for the fitting parameter c(L,q). 



L 


q = l 


q = 2 


q = 3 


q = i 


32 


4.76(3) 


4.29(4) 


3.97(3) 


4.08(4) 


48 


4.75(3) 


4.25(1) 


4.00(2) 


4.06(2) 


64 


4.70(2) 


4.23(2) 


4.00(3) 


4.09(2) 


80 


4.71(6) 


4.19(3) 


3.99(3) 


4.12(2) 


128 


4.54(8) 


4.26(4) 


3.96(3) 


4.17(2) 



TABLE XI: Results of the approximation for the fitting parameter d(L,q). 



L 


q=l 


q = 2 


q = 3 


q = A 


32 


-0.07(4) 


-0.30(7) 


-0.29(5) 


-0.26(7) 


48 


-0.10(4) 


-0.21(3) 


-0.35(3) 


-0.14(4) 


64 


-0.05(3) 


-0.19(4) 


-0.34(6) 


-0.13(3) 


80 


-0.06(9) 


-0.11(5) 


-0.33(6) 


-0.13(3) 


128 


0.14(12) 


-0.23(7) 


-0.31(6) 


-0.08(3) 



VI. DISCUSSION 



Using the dual lattice (see Fig. |2J) allows us to avoid finite size shift of the critical point 
for the bond percolation and to diminish it for spin models. The accuracy of definition of 
the critical point on the finite lattice play the principal role for the investigation of the tails 
scaling. The high quality of our approximation is proved by remarkable symmetry of the 
crossing probability with respect to the critical point p c . In Fig. EJand Fig. [7|we can observe 
that the two branches p — p c > and p — p c < practically coincide. 

The two different scaling region of the crossing probability clearly seen in Fig. ^i) can 
explain the long time uncertainty about its shape. In Ref.jljJ the scaling index for the 
percolation threshold for 2d percolation model was found £ ~ 1.9(1). This result coincides 
with our approximation of the body of the crossing probability for percolation ( ~ 1.864(12). 



In more recent works 



ll 16, 12, H, 



21] the tails region for percolation was investigated 



which is described by the scaling formula _Dexp(— cL(p — p c ) u 4 / 3 ). The crossover form to 
Gaussian-like behavior to slope 4/3 is observed in figures of Ref . [22I] . It seems, near the 
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ated 



critical point the behavior of the crossing probability is parabolic. The rounding happens 
in the interval r < 0.1. This interval is relatively small in comparison with regions of the 
body 0.1 < \t\ < 1.5 and tails 1.5 < |r| < 4 as it can be seen in Fig. 

We have five fitting parameters a(q), b(q), z(q), c(q), d(q) in expressions (fTU|) and (jlljl . In 
Fig. EK)-Fig. 01) we see the crossover region between the body and the tails of the crossing 
probability. In this region the function f 3 touched the line f±. It means, that in some point 
To the values of the functions are equal, therefore 

- a(L, q) + b{q)4 {q) = -d(L, q) + c(q)r , (15) 

and first derivatives of this functions are equal too 

z{q)b{qK {q) - 1 = c(q). (16) 

Substituting c(q) from Eq. (JTBJ) in Eq. (fTH|) we obtain expression for b(q) 

d(q) - a(q) 

(*(?) - iK ( 

If the crossing probabilities at the critical points A(q) (and logarithms a(q)) can be calcul 
analytically by conformal field methods (at least for the percolation it is possible) 

Baa 

then only four independent parameters b(q), T (q) and z(q), a(q) remain for the crossing 
probability. 

The main statements for the crossing probability 717^ are: 

• In accordance with scaling theory the finite size scaling of the crossing probabilities 
may be eliminated by introducing the scaling variable r = L(p — p c (L, q)) v<yq \ The 
crossing probability as a function of r does not depend on the lattice size L. 

• The body of the crossing probability scales vr/j S (r) ~ A(q) exp(— b(q)r z( - q ^). 

• The tails of the crossing probabilities scales iihsij) — D(q) exp(— c(g)r). The value of 
parameter D(q) is about 1. 

• The finite size scaling for q = 4 does not describe by the analytical value of the 
correlation length index v(q = 4) = 2/3. We obtain some scaling index x(q = 4). This 
index x(q = 4) ~ 0.759(4) for the tails of the crossing probability (see Table |^) or 

~ 0.728 for the body of the crossing probability (see Table ITT]) . 
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